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Abstract 

This paper introduces a matrix analog of the Bessel processes, taking values in 
the set E of real square matrices with nonnegative determinant. They are related to 
the well-known Wishart processes in a simple way: the latter are obtained from the 
former via the map x i— > x T x. The main focus is on existence and uniqueness, questions 
which lead us to develop several new results of potential theoretic nature concerning 
the space of real square matrices. Specifically, the function w(x) = | dct x\ a is a weight 
function in the Muckenhoupt A p class for — 1 < a < (p = 1) and — 1 < a < p — 1 
(p > 1). The set of matrices of co-rank at least two has zero capacity with respect 
to the measure m{dx) = \ detx\ a dx if a > —1, and if a > 1 this even holds for the 
set of all singular matrices. As a consequence we obtain density results for Sobolev 
spaces over E° with Neumann boundary conditions. The proofs of these results rely 
crucially on the stratification of the state space into smooth manifolds consisting of 
fixed-rank matrices. 
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1 Introduction 

The Wishart processes constitute a class of matrix- valued generalizations of the squared 
Bessel (BESQ) processes, and were studied systemically for the first time by Bru [I]. 
Her work has subsequently been extended in various directions, for instance by Donati- 
Martin et al. [7] and Cuchiero et al. [6], among others. The fact that the BESQ processes 
have well-behaved matrix-valued analogs begs the question of whether the same is true 
for the Bessel (BES) processes. Since Wishart processes take values in S+, the cone 
of positive semidefinite matrices, a natural attempt would be to consider their positive 
semidefinite square roots. However, as was already remarked by Bru [1], these processes 
appear intractable. Our goal is to show that by passing to a larger state space, a more 
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well-behaved class of processes can be obtained. Specifically, the state space will be the 
set 

E = {x G M d : detx > 0}, 

where M rf is the Euclidean space of all d x d real matrices, endowed with the usual inner 
product x • y = tv(x T y) and norm ||x|| = ^/x~»y. The matrix-valued Bessel process with 
parameter 5 > and matrix dimension d, abbreviated BESM(J, d), will be an E- valued 
Markov process whose generator (acting on a suitable function space) is given by 

A = ^A+ 6 -^x- T .V, (1) 

where x~ T = (x~ 1 ) T , V is the d x d matrix with elements d Xi . (so that for / G C 1 (M d ), 
V/ is its gradient), and A = V • V = ^-x y is the Laplacian. Notice that A is the 
generator of the BES(£) process when d = 1. A general discussion on BES and BESQ 
processes is available in [191 Chapter XI], and a specialized treatment of the BES(<5) 
process, < 5 < 1, can be found in [SJ. 

The singular behavior of the drift part of A makes a treatment via the associated martin- 
gale problem inconvenient. Instead we use the machinery of Dirichlet forms, which turns 
out to be very suitable for proving existence. We will apply the theory as presented in the 
book by Fukushima et al. [12]. The crucial fact here is that A is a symmetric operator with 
respect the measure m(dx) = | det x\ s_1 dx (where dx is Lebesgue measure on Nl d ), which 
is a consequence of an integration by parts formula (Theorem [2]). The precise definition 
of the BESM process and semigroup, together with the discussion of existence, is given in 
Section [2J 

Uniqueness, discussed in Section HI is a much more delicate issue, whose resolution oc- 
cupies the majority of the paper. We will obtain Markov uniqueness of the associated 
semigroup, relying on density results for Sobolev spaces with Neumann boundary condi- 
tion (Theorem [8]) . In the process we are lead to prove several results, interesting in their 
own right, about the measure m and its interaction with the state space. Specifically, we 
prove that | detx| a is locally Lebesgue integrable on M rf precisely when a > — 1 (Theo- 
rem [1]), and that it is a weight function in the Muckenhoupt A p class when a G (—1,0] 
and p = 1, and when a G (— l,p — 1) and p > 1 (Theorem [9} . This exactly parallels the 
situation for the weight function on R given by t a . Moreover, we show that the matrices 
of co-rank at least two form a set of zero capacity with respect to m, and that the set of 
all singular matrices has zero capacity if 5 > 2 (Theorem I10p . The latter result relies on 
the fact that the sets 

M k = {x G M d : rank x = k} 

are smooth manifolds forming a stratification of NL d . In particular, it allows us to deal 
with the highly non-convex, non-Lipschitz structure of the state space E. In fact, the 
interior E° does not even lie on one side of its boundary dE, as can be seen by considering 
the lines {tx : t G R}, x G M d (which intersect dE): If d is even, each line lies either 
entirely inside E, or entirely outside E°. In the future we expect similar techniques to 
be useful for analyzing a wide class of processes on semialgebraic (or, more generally, 
stratifiable) state spaces. Some of the groundwork for this is laid in [S]. 
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We emphasize that the case < 5 < 2 is the difficult one, due to the fact that the boundary 
of the state space is reached almost surely. This is analogous to the scalar Bessel processes. 
When 5 > 2, on the other hand, the boundary of E is never reached, and the BESM process 
can be realized pathwise as the strong solution to a suitable stochastic differential equation. 
This is discussed at the end of Section [3l 

Let us say something about why the BESM processes are natural analogs of the BES 
processes, other than the resemblance of their infinitesimal generators and the fact that 
they coincide for d = 1. The main reason is that the process X T X, where X is BESM(<5, d), 
is a Wishart process with parameter a = d — 1 + 5, denoted WIS (a, d). As a consequence, 
||X|| is BES (da), and det X is a time-changed BES (5) process. These results are presented 
in Section It is interesting to note that the restriction a > d — 1, which corresponds 
exactly to the range where the symmetrizing measure m has the Radon property, is where 
the stochastic differential equation for the Wishart process has a weak solution, see [H 
Theorem 2]. 

A second motivation, which was the original "clue" that led us to consider the generator 
A, is as follows. Let X be an M d -valued Brownian motion starting from I (the identity 
matrix), and let tq be the first time det At hits zero. Then det Xt/\ TQ is a martingale, and 
we may use it to change the probability measure. An application of Girsanov's theorem, 
using the identity Vlndet(x) = x~ T , shows that under the new measure, the process 



is Brownian motion. Thus X becomes a Markov process whose generator is A with 
(5 = 3, and its determinant is positive by construction. This is fully analogous to the 
well-known construction of the BES (3) process as Brownian motion "conditioned to stay 
positive" . 

In addition to the notation already introduced above, the following conventions will be in 
force throughout the paper. 

• For any subset U C M rf , C(U) denotes the space of real- valued continuous function 
on U, equipped with the relative topology. The subspace of compactly supported 
functions is denoted by C C (U). If the functions take values in a space V, we write 
C(U; V). A function on U is called k times continuously differentiable (k G NU{oo}) 
if around each x G U it coincides with some k times continuously differentiable 
function on M rf . The space of such functions is denoted C k (U). The meaning of 

(£/), C^(U; V), etc. should be clear. The following remark is trivial but important: 
the notion of "compactly supported" and "vanishing at infinity" depends on the class 
of compact sets, which in turn depends on U. In particular, if U is open in M rf , a 
set K C U is compact in U if and only if it is bounded, closed in M rf , and bounded 
away from dU. If on the other hand U is closed in M d (as is the case when U = E), 
compact subsets need not be bounded away from dU. The same remark applies to 
any property that holds locally on U. 

• As already mentioned, the sets Mk, k = 1, . . . ,d, constitute a stratification of M rf 
into smooth manifolds, and the dimension of M\. is d 2 — {d — k) 2 . Note that dE = 




t > 
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Ufc<d_iMfc. The tangent space of at x S is well-defined, and for any v £ M 
we say that v is tangent to dE at x £ if v lies in the tangent space at x of 
the manifold Mk that contains x. We refer to |18| for background on differential 
geometry. 

• 0(d) is the orthogonal group over R d , and T(d) is the group of upper-triangular 
d x d real matrices with strictly positive diagonal entries. The set of nonsingular 
matrices (i.e., the general linear group) is homeomorphic to 0{d) x T(d) via the 
QR-decomposition. 

• For x 6 M d , we let adj x denote the adjugate matrix of x (i.e., the transpose of 
the matrix of cofactors). We will freely use the identities xadjx = (detx)I and 
Vdet(x) = adj x, as well as the fact that adj x = if and only if rankx < d — 2. 



2 Definition and existence 

The definition of the BESM process departs from the differential operator A in ([I]) acting 
on functions in C? j_ , where 

C l x = {/ e C c( E ) : X ~ T • W is bounded} . 

As we will see momentarily, (A, C^j_) is a symmetric operator on L 2 (E,m), where the 
measure m is given by 

m(dx) = | det x\ s ~ 1 dx. 

Occasionally m will be viewed as a measure on M rf , or on subsets other than E. It is 
crucial that m is a Radon measure on M d , and hence on i?, when 5 > 0. The proof of this 
fact relies on the following lemma, which is an elegant consequence of the group structure 
associated with the QR-decomposition of an invertible matrix. 

Lemma 1 Let f : M d — > R be nonnegative and measurable. Then 

[ f(x)\detx\- d dx= I f(QR)n(dQ)HR^dR, 

JM d JO(d)xT(d) ~~[ 

where fi is a non-normalized Haar measure on 0(d), and dR = Ili<j dRij. 

Proof. If / (x)\ det x\~ d is integrable, the result is exactly Proposition 5.3.2]. The 
general case follows by monotone convergence. ■ 

The Radon property of m can now be proved. 

Theorem 1 Let a € R and define w(x) = \ detx| a . The function w is locally integrable 
on M d if and only if a > — 1 . 

Proof. Let A C M d be relatively compact. Since dE is a nullset, we may assume that 
An dE = 0. Then there is a cube K C T(d), say K = Y\ iK Lij, of bounded open intervals 
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Iij such that la C (0, oo) and Ij,- C R (i < j), and such that A C 0(d) • K. Hence, the 
change-of- variable formula in Lemma [1] yields 

j \ det x\ a dx < n(0{d))^Y[J^ dR^i^J | det R\ a+d JJ R^dR, 

=Ko(d))(n[ dR t Mn[ R^ d - i dR t 



d 

i<j " 1 ij ' V i = l T 

where fi is a non-normalized Haar measure on 0(d). The right side is finite, provided 
a > —1. If on the other hand a < — 1, take 1^ = (0, 1) for all i < j, and set A = 0(d) ■ K. 
Then A is relatively compact, but f A w(x)dx = oo. ■ 

Consider the following differential operator: 

V*G = -(V + (5- l)x~ T )»G, G G C 1 (E;M d ). 

Notice that A = — ^V*V. With this in mind, the integration by parts formula below 
shows that A is indeed a symmetric operator — indeed, taking G = Vg with g G ± , it 
implies that the equality 

(/, -Ag) = l -(Vf,Vg) 

holds for all / G C^(E). Here and in the sequel, (•,•) denotes the inner product on 
L 2 (E,m). Specifically, we write 

= / f{x)g(x)m(dx) and (F,G}= / F(x) • G(x)m(dx), 
Je Je 

where /, g G L 2 (E,m) and F,G G L 2 (E,m;M. d ). The overlapping notation should not 
cause any confusion. 

Theorem 2 (Integration by parts formula) Assume that 5 > 0, and let f G Cl(E) 
and G G C 1 (E; M d ). If 5 < 1, assume that G(x) is tangent to dE at x for all x G dE. If 
5 < 1, assume in addition that G(x) • x~ T is locally bounded. Then 

(V/,G) = (/,-V*G). 

Proof. See Section [5j ■ 

Remark 3 It is not hard to show that local boundedness of G(x) • x~ T implies that G(x) 
is tangent to dE at x for all x G dE. 

The BESM process can now be defined as follows. 

Definition 1 (BESM semigroup) A strongly continuous sub-Markovian contraction semi- 
group (T t :t> 0) on L 2 (E,m) is called a BESM(5,d) semigroup if its generator extends 

iAC 2 ,). 
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An E- valued Markov process X is called m- symmetric if its transition function pt(x,dy) 
is m-symmetric. In this case the operators f E f(y)pt( x ,dy), where / is bounded and in 
L 2 (E,m), can be extended to all of L 2 (E,m), see JT2J page 30]. This extension is called 
the L 2 (E,m) semigroup of X. 

Definition 2 (BESM process) An E -valued m-symmetric Markov process whose L 2 (E,m) 
semigroup is a BESM(5, d) semigroup is called a BESM(5, d) process. 

A substantial portion of this paper is devoted to the question of uniqueness of the BESM 
semigroup and process. Existence, on the other hand, is straightforward via the theory 
of Dirichlet forms. In view of the discussion so far it is natural to introduce the bilinear 
symmetric form 

£(/,</) = ^<V/,V 5 ), f,geCl(E). 

This form is closable on L 2 (E, m), as can be deduced from Theorem[2]as follows. Pick a se- 
quence (/„) in Cl(E) converging to zero in L 2 (E, m), such that lim m>n £(f n -f m , f n -fm) = 
0. We must show lim n £(f n , f n ) = 0. Since (V/ n ) is a Cauchy sequence in L 2 (E,m;M d ), it 
has a strong limit F. For any $> G C£°(E; M rf ) vanishing near dE, Theorem [2] yields 

{F *> = lim (V/ n , $) = lim(/ n , V*$) = 0. 

n— >oo n 

It follows that F = m-a.e., as required. We denote the closure of (£, C^(E)) by 

(£,D(£)). 

An application of [12^ Theorem 3.1.2] then shows (after routine verification of the condi- 
tions of that theorem) that (£,D(£)) is a Dirichlet form. Moreover, it is 

(i) regular: D(£)nC c (E) is || * ||oo~dense in C c ( y E') and f^-dense in Z)(£), where 

£i(f,f) = (f,f) + £(fj). 

(ii) strongly local: for every f,g& D{£) with compact support such that g is constant 
on a neighborhood of the support of /, we have £(f,g) = 0. 

The integration by parts formula implies that the equality £(f,g) = {—Af,g) holds for 
all / G C c 2 j_ (indeed, for all / G C 2 if 5 > 1) and all g G Cl{E), hence for all g G D{£) by 
continuity. It follows that the generator associated with £ coincides with A when acting 
on functions / as above. With some abuse of notation, we therefore let 

(A,D(A)) 

denote the generator of £, noting that the domain D(A) contains C 2 j_, and even contains 
C 2 (E) if 5 > 1. In particular, the semigroup (T t : t > 0) on L 2 (E,m) associated with £ 
and A is a BESM(<5, d) semigroup. 

A corresponding BESM(<5, d) process is then obtained as the m-symmetric Hunt process 
X on E corresponding to the Dirichlet form {£ , D(£)), see |12( Theorems 7.2.1 and 7.2.2]. 
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The strongly local property of £ implies that this process has continuous paths. However, 
it is not guaranteed a priori that X is conservative; we now prove that it is, thereby 
obtaining existence of the BESM(<5, d) process. In the following, let ¥ x be the law of X 
starting from x G E. 

Proposition 1 The Dirichlet form £ is conservative, i.e. the semigroup (Tt : t > 0) 
satisfies T^l = 1 for all t > 0. Consequently, X can be chosen so that 

F x (X t G E for all t > 0) = 1 for all x G E. 

Remark 4 For / G L°°(E,m) with / > 0, one defines Ttf = \im n Ttf n , where f n G 
L 2 (E,m) and / n t / wt-a.e. This definition is independent of the choice of sequence (f n )- 

Proof. By [121 Theorem 1.6.6], £ is conservative if there is a sequence (/ n ) C D(£) such 
that < f n < 1 and lim n f n = 1 m-a.e., and such that 

lim holds for any g G n L l (E,m). 

n 

To construct such a sequence, let G C 2 (R) satisfy 0(f) = 1 for t < 0, = for t > 1, 
and |<^'(i)| and |</>"(i)| uniformly bounded. For n > 1 define 



/ n (x) = - n) 



Differentiating twice yields 



.r 



V/„(a;) = ^(||x||-n) 

I 



|x|| 

Afn{x) = <f>"(\\x\\ - n) + 0'(||x|| - n)- 

In particular we have / n G C? j_ C -D(»4). Since also V/ n and A/ n both vanish outside 
the set E n = {x G E : n < \\x\\ < n + 1}, we obtain 



|£(/n,3)l 



Af n (x)g(x)m(dx) 



<esssup xeEn \Af n (x)\ / |g(x)|m(dx), 



where Holder's inequality was applied. For n > 1, the essential supremum is bounded by 
a constant c > that is independent of n. Hence 



^2\£(fn,9)\<cJ2 [ \9 

n>l n>0 En 



(x)\m(dx) = c\\g\\ L i {E)m) < oo. 



We deduce that lim n \£(f n ,d)\ = 0, showing that £ is conservative. The statement about 
X now follows from [121 Exercise 4.5.1]. ■ 

We finish this section with a preliminary result regarding C 2 j_ . In Section UJ this set will 
be studied further. 

Proposition 2 The set C 2 j_ is a core for £. 

Proof. We apply a slightly modified version of (SJ Theorem 5.1], where in the proof we 
use Theorem 6.2 instead of Theorem 2.3 (both in the same reference.) This shows that 
C? j_ is dense in C^(E) with respect to both || • and £\. ■ 
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3 Relation to the Wishart process 



We now describe the properties of the transformed processes XX, \\X\\, detX, where X 
is a BESM(<5, d) process realized as an m-symmetric Hunt process with BESM(5, d) semi- 
group, as discussed in the previous section. We always take 5 > 0. The main observation 
is the following. Let 

<f> :M d S d + , x^x T x, 

and define 

A wls g{z) = tr(2zV 2 + (d - 1 + S)V)g(z), z G S d + , g G C C °°(S^). 

This is the generator of the WIS(a, d) process with a = d — 1 + 8, see [3]. Now, for any 
g G C C °°(S^) and any x G E, one readily verifies the identities 

V(g o $)(x) = 2xVg(x T x), -A(g o $)(a;) = tr(2x T xV 2 + dSJ)g{x T x). 

Consequently we have 

go$eD(A) and A(g o $) = {A W1S g) o $ G C 6 °°(S^). (2) 

This relation will let us deduce that X T X is a Wishart process. The proof uses the 
following standard result. 

Lemma 2 Pick any f G D(A) such that Af is locally m-integrable on E. For quasi-every 
x G E, the process 

f(X t )-f(x)- [ Af(X s )ds, t>0, (3) 
J o 

is a square integrable martingale under F x , and Jq \ Af(X s )\ds < oo for all t > 0, ¥ x -a.s. 

Proof. This can be deduced from |12|, Theorem 5.2.2, Corollary 5.4.1, and Theorem 5.1.3]. 
The local integrability assumption is to ensure that \Af{x)\m(dx) is a Radon measure, 
which is needed in |12[ Corollary 5.4.1]. ■ 

Remark 5 The exceptional set for which the conclusion of the lemma fails depends on 
the function / in general. However, while it will not be proved it here, we conjecture that 
X can be chosen so that its transition function is absolutely continuous with respect to m. 
This implies in particular that if / and g are measurable and satisfy / = g, m-a.e., then 
the equality f(X t ) = g(X t ) holds F x -&.s. for every x G E and every t > 0. By a standard 
argument using |10l Proposition 1.1.5], where the Banach space L is taken to be L 2 {E, m), 
one can then show that ([3]) defines a P x ' martingale for every / G D(A) and every x G E 
such that E x [/q \ Af{X s )\ds\ < oo, t > 0. In particular the "quasi-every" caveat can be 
removed in this case. 

Theorem 6 Let X be a BESM(<5, d) process as above. 

(i) X T X is a WIS (a, d) process, where a = d — 1 + 5. 

(ii) \\X\\ is a BES(a) process, where a = d(d — 1 + 5). 
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(Hi) detX is a time-changed BES(5) process. More specifically, define 

a(t) = / || adj X s \\ 2 ds, £ u = detX 6(u ), 
J o 

where b(u) = inf{i > : a(t) > u} is the left- continuous inverse of a. Then £ is a 
BES(5) process stopped at a(oo). 

Remark 7 The precise meaning of part (i) of the theorem is that for all x outside some 
exceptional set N, the law of X T X under ¥ x is that of a WIS (a, d) process. The statements 
(ii) and {Hi) should be interpreted analogously. Due to Remark [5l absolute continuity of 
the transition function would let us take N to be empty. 

Proof. Lemma [2] together with ([2]) imply that for any countable family C C C^°(S^) 
there is an exceptional set N C E such that ([3]) is a square integrable martingale under 
F x for all x <£ N and all / = g o $ with g G C. Setting = Xj X u it follows that the 
process 

g(Z t ) - g(x T x) - [ A wls g(Z s )dx, t > 0, 
Jo 

is a martingale under P x for all g G C Using standard arguments (see for instance [20^ 
Theorem V.20.1]) this lets us deduce that X T X is a weak solution to the stochastic 
differential equation 

dZ t = y/ZtdWt + dWj sfz t + (d - 1 + 6)ldt, Z = x T x, 

for all x ^ N, where W is M rf - valued Brownian motion under P x . (This also uses the 
well-known correspondence between the above stochastic differential equation and „4 W .) 
Hence (i) follows. Part (ii) is immediate from the well-known fact that the trace of a 
WIS(a, d) process is a BESQ(ad) process. For Part (Hi) we set Z = X T X and note 
(see [U Section 4] ) that det Z satisfies 

det Zt = det Zq + 2 \J det Z s tr(adj Z s )d(3 s +81 tr(adj Z s )ds 
Jo Jo 

for some standard Brownian motion /3. Defining 

a- x (u) 



Pu = / ^tr(ad} Z s )d(i s , 



which is Brownian motion on [0, a(oo)), and using the identity adj z = (adjx) T adjx, we 
get 

det Zt = det Zq + 2 / \j det ~Z~ s dfi a t s \ + 5a(t). 

Applying the time-change it = a(t) it follows that det Zu.\ satisfies the stochastic differen- 
tial equation for the BESQ(5) process, stopped at a(oo) since b is constant on [a(oo),oo). 
Since detXt = y/ det Zt , the result follows. ■ 

By a similar argument one can show that X itself is a weak solution to a certain stochastic 
differential equation when 5 > 1. Indeed, we saw in Section [2] that C^(E) C D(«4) if 5 > 1, 
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and Theorem [I] implies that 1/ det (a;) is locally m-integrable in this case. We then deduce 
from Lemma [2] that ([3]) is a F x martingale for all / G C C C^(E), C countable, and all x 
outside some exceptional set N. This in turn implies that for all x ^ N, X satisfies 

X t = x + W t + ^1 [ X~ T ds, t > 0, 
* Jo 

where W is Brownian motion under ¥ x . In particular, X is a semimartingale. On the 
other hand, if < 5 < 1, X is not a semimartingale — if it were, so would detX, which 
however is a time-changed BES(<5) process by Theorem [6] and hence not a semimartingale. 
In the case 5 = 1 one can show that X is a semimartingale satisfying a Skorohod-type 
stochastic differential equation. We do not elaborate on the details here. 

We close this section with a pathwise construction of the BESM process as the strong 
solution to a stochastic differential equation. This construction only works when 5 > 2 
and the process starts from E° . We conjecture that the following result remains true for 
all 5 > 1. 

Proposition 3 Suppose 5 > 2, and let W be standard NL d -valued Brownian motion de- 
fined on some probability space. The stochastic differential equation with state space E° , 

dX t = dW t + ^X; T ds, (4) 

has a unique strong solution. 

Proof. Since x i— > x~ T is locally Lipschitz continuous on E°, a strong solution X exists 
on [0,0, where £ = lim n inf{i > : det Xt < n" 1 or ||X^|| > n}. We claim that £ = oo. 
To see this, set Z = X T X and note that we have 

Z t = f X s dW s + f dWjxJ + (d - 1 + S)lt, t < C, 
Jo Jo 

where we used the equality dW t dW t = dldt. Defining W t = f*(X] 'X,)" 1 / 2 X s dW s , t < £, 
we have t 

Z t = I ^/Y s dW s + [ dWjy/z~ s + (d - 1 + 8)It, t < c, 
Jo Jo 

and after verifying that W is again M^-valued Brownian motion on [0,0, it follows that 
Z is a WIS(d — 1 + S,d) process on [0,0 (this calculation is of course closely related to 
the one in the beginning of Section [3j) Since 5 > 2, this implies that detXt = \J det Zt 
stays strictly positive, and that ||^t|| 2 = tv Zt is nonexplosive. Hence £ = oo as claimed, 
and the result is proved. ■ 

4 Uniqueness 

The goal of this section is to establish uniqueness of the BESM semigroup. Specifically, 
we will prove that (A, C 2 j_) is Markov unique. This means that there is at most one (and 



10 



hence exactly one) symmetric sub-Markovian strongly continuous contraction semigroup 
on L 2 (E,m) whose generator extends (A,C 2 see [HI Definition 1.1.2]. Since the Hunt 
process corresponding to such a semigroup is unique up to equivalence^, this form of 
uniqueness will hold for any realization of the BESM process as a Hunt process. 

Note that m(dE) = 0. Therefore L 2 (E,m) and L 2 (E°,m) can be identified, implying that 
it is enough to prove Markov uniqueness of (A, C 2 ± ) as an operator on the latter space. 
To do this we will apply a general result by Eberle [9l Corollary 3.2] that relies on studying 
the relationship between various weighted Sobolev spaces, which we now introduce. To 
simplify notation we henceforth write 

= E° = {x G M d : det x > 0}. 

Observe that l/(det x) s ~ 1 is locally integrable on Q. Hence by \17\ Theorem 1.5], L 2 (fi, m) 
is continuously imbedded in L 2 (Q,dx). In particular, every / G L 2 (Q,m) has a unique 
distributional gradient Df, and one can define the weak Sobolev space 

W l ' 2 {n,m) = {/ € L 2 (n,m) : Df G L 2 (Q, m; M d )} . 

Equipped with the norm 

1 /2 

ll/ll^(n,m) = (J \f(x)\ 2 ™(dx) + J \\Df(x)\\ 2 m(dx)^ , 

W 1,2 (£l,m) becomes a Hilbert space \17\ Theorem 1.11]. We also consider the following 
strong Sobolev spaces (here the word completion is always meant with respect to the norm 

II " llvi /1 ' 2 (n,m)) : 

H^^m) = completion of C°°(ft) 

£Tj_' 2 (r2, m) = completion of C 2 j_ 

#o' 2 (n,m) = completion of C c °°(^). 

These are all Hilbert spaces by construction, and we automatically have 

# 1,2 (ft, m) C H 1 /^, m) C F 1,2 (fi, m) C W 1 ' 2 ^, m). (5) 

In the present section we will expend significant effort to show that for any 5 > 0, the 
last two inclusions are equalities; and if <5 > 2, all three inclusions are equalities. This 
will lead to Markov uniqueness of the BESM semigroup. Specifically, we will prove the 
following: 

Theorem 8 The following statements hold, 
(i) IfO < 6 < 2, then F^ 2 (0,m) = W l > 2 {Q,m). 

1 Two symmetric Hunt processes are called equivalent if their transition functions coincide outside a 
properly exceptional set, see Section 4.1 in [12] . 
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(ii) IfS>2, then H^ 2 (n,m) = W 1 ' 2 (9,,m). 
In either case, (A,C 2 ,) is Markov unique. 

The proof of this theorem requires some preparation. As in the one-dimensional situation, 
the case < 5 < 2 is the most delicate. The following result, interesting in its own right, 
is crucial for handling this case. We let \A\ = f A dx denote the Lebesgue measure of a 
measurable subset A C M^. 

Theorem 9 (Muckenhoupt property) Let a G R and define w(x) = | detx| a . 

(i) If — 1 < a < 0, then w lies in the Muckenhoupt A\ class. That is, there is a constant 
C > depending only on d and a, such that 



for every ball B C M d . 

(ii) If —1 < a < p — 1, p > 1, then w lies in the Muckenhoupt A p class. That is, there 
is a constant C > depending only on d, a and p, such that 



for every ball B C M d . 

Proof. The proof of part (i) is given in Section [SJ Part (ii) follows from (i) via [21\ 
Proposition IX. 4. 3]. ■ 

An application of this theorem with a = 5 — 1 shows that | detx) 5-1 is an A2- weight if 
< 5 < 2. In particular, [161 Theorem 2.5] then shows that the last inclusion in ([5]) is in 
fact an equality. 

The other key ingredient in the proof of Theorem [8] is to show that the sets of low-rank 
matrices are small. To measure the size of a set we use the following notion of capacity 
(referred to as the (l,m)-Sobolev capacity in [13]). For any subset A C M rf , define 



where the infimum is taken over all / G H 1 ' 2 (M d , m) with / = 1 on an open neighborhood 
of A, see \13\ Definition 2.35]. We have the following result: 

Theorem 10 Let 5 > 0. For k G {0, . . . , d - 2}, we have Cap(M fc ) = 0. If 5 > 2, the 
same thing also holds for k = d — 1. 

The core of the proof of Theorem [10] is an application of the following lemma, which 
bounds the growth of the determinant function near a point x G M^. 

Lemma 3 Let k < d — 1. There is a locally Lipschitz function c/% : — > R + such that 






det(x + v)\ < Ck(x)\\v 



d-k 



x G M k , v G M' 



d 



v\\ < 1. 



12 



Proof. By [3j Corollary 5], 



d 

|det(x + 4 -det(x)| < ^p^a^x), . . . , a d (x))\\v\\ % , 

i=l 

where cr(x) = (cri(x), . . . ,ad(x)) is the vector of singular values of x, and pi is the i:th 
elementary symmetric polynomial in d variables. Now, pd-i{cr\(x) , . . . , a^(x)) consists of a 
sum of terms, each of which is the product of d—i distinct elements of o~(x). However, since 
ranks; = k, only k of those elements are nonzero. Therefore the product must contain 
at least one zero factor whenever d — i > k, implying that pd-i{o-\{x), . . . ,<Jd(x)) = for 
these i. Since in addition \\v\\ < 1 and detx = 0, we get 

d 

\ det(x + v)\ < \\v\\ d - k Pd-i((^i(x),...,a d (x)). 

i=d—k 

The local Lipschitz property follows from the smoothness of pd-i and the fact that the 
singular value map is Lipschitz continuous, see \\.4\ Theorem 7.4.51]. ■ 

In proving Theorem [TUJ the case k = d — 1,5 = 2 turns out to require separate treatment 
using the following lemma. 

Lemma 4 For each e < 1 there is a Lipschitz function <j) £ : R+ — > R such that < <p e < 1, 
4> £ = on [e, oo), (j) £ = 1 on a neighborhood of zero, and 

lim/ \<f>' £ (t)\ 2 tdt = 0. (6) 
e l° Jr + 



Proof. Define functions g e and h £ on R + by 



9e(t) 



and h £ (t) = < 



' {t/ef t 6 [0,e 1+1 / £ ) 

2e-e- x Ih t e [e l+l / £ ,2e l+l,£ ) 
t 6 [2e 1+1//£ , oo) 



We claim that the function 



g £ + h £ has the stated properties. It is not hard to check 



that < <f> £ < 1 and that <ft £ equals zero on [e, oo) and one on [0,e 1+1 / £ ). The Lipschitz 
property then follows easily. It remains to verify ([6]). First, note that 



/ 



\g' £ {t)\ 2 tdt = e 



2-2e 



dt = -^o (eio). 



Moreover, since \h' £ \ = \g' £ \ on [0, e 1+1 / £ ), and since 

\h' £ (t)\ 2 tdt = e-^^ 1—£ >- = &^0 (40), 

;l+X/e ^ ^ 

it follows that lim £ | Jr \h' £ {t)\ 2 t dt = 0. We now deduce (jSJ). ■ 
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The proof of Theorem [10] exploits the fact that is a smooth manifold — specifically, we 
need the notion of an open tube segment. To introduce this notion, let M be any smooth 
ni-dimensional embedded submanifold of R n (ni < n), and fix any point x € M. Then 
take an open neighborhood U C R n of x, and a diffeomorphism 

$ : U -> $(£/) 

such that $(Mn*7) = M x {0} G R™ 1 x R™ 2 , n 2 = n-ni, for some open subset M C R ni 
containing $>(x) (these objects exist by definition of smooth manifold). Let G be an open 
ball in R ni , centered at with its closure contained in M. Let B £ be the open ball 

in R n2 , centered at the origin, with radius e. For all sufficiently small e, say e < £0) the 
closure oi G x B e lies in &(U). Define the set 

U £ = <& _1 (G x B e ). 

We refer to t7 £ as an open iu6e segment around M centered at x. It follows from the 
above that for each e < £q, U £ is an open subset of R n containing x {0}). In our 

setting, R n is identified with M rf , so we have n = d 2 , and M is identified with for the 
particular k under consideration, so that r&2 = (d — k) 2 and n\ = d 2 — (d — k) 2 . 

We now have all the ingredients needed for the proof of Theorem [TU1 

Proof of Theorem 1101 Let U e be an open tube segment around centered at some 
fixed x, as described above. Let ir : R ni x R" 2 — > R™ 2 be the coordinate projection onto 
the last 7i2 coordinates. Let if £ C°°([0, oo)) be a cutoff function valued in [0, 1], equal to 
one on [0, 1/2], equal to zero on [l,oo), and with |y'(£)| < 3 for all t > 0. Define 

g £ (x) = if {e- 1 \\tt(x)\\) , xeGxB eo , 

and then 

g e (x) =g e o$(x), x£U £0 . (7) 

It is clear that g £ lies in C°°(U £0 ), vanishes outside U e , and is equal to one on an open 
neighborhood U £ of $ _1 (G x {0}). Let us compute Vg e . Pick any v G and any 
x £ U £o . By the chain rule, 

Vg £ (x) mv = Vflf e ($(x)) • &(x)[v], 

where $'(x)[f] denotes the Frechet derivative of <& at x applied to v. The chain rule also 
gives, for x G G x i? £o , 

Vff £ (J) .v = h(f (e^llTrCxJII) • tt'(x)H. 

The derivative of tt is given by 7r'(x)[f] = ir(v), and since orthogonal projections are 
contractions, we deduce 

l|Vflfe(x)|| < -, 
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and consequently 

l|V 5e (x)|| < -C, 

where C = sup a;gt / eo ||$'(a;)|| (operator norm) is finite due to the fact that the (compact) 
closure of U £0 is contained in a set where 3> is smooth. It follows that 

2- 



Ue 



(\g £ (x)\ 2 + \\Vg £ (x)\\ 2 ) m(dx) < (l + m(U £ ). (8) 



Suppose we can show that the right side tends to zero as e 1 0. Then, multiplying g £ by 
a smooth function ip that is zero outside U £0 and one on some open neighborhood U' e 
of x, we obtain functions f £ G .ff 1,2 (M , m) certifying that Cap(Mfc n U £Q ) = 0. Since 
Mj, is covered by a countable union of such sets U' e , |13l Theorem 2.37] implies that 
Cap(Mfc) = 0, as required. 

Thus, it remains to show that e~ 2 m{U £ ) — > as e — > 0. A change of variables yields 
m(U £ )= / \detx\ s ~ 1 dx= / | det v)) 5 " 1 J(y, v)dy ® dv, 

where J(y, v) is the Jacobian determinant, and where we write R n = R ni x R™ 2 ,x = (y,v), 
dx = dy <8> dv (Lebesgue measures). Since <I> -1 is smooth on a set containing the closure 
of G x B £ , it is Lipschitz on G x B £ with some constant K\. Moreover, J(y,v) is bounded 
by some constant k 2 . Together with Lemma [3] (and the fact that <3? _1 (y,0) G M k ), this 
yields 

m(U E )<K 2 [ \det($- 1 (y,0) + <l>- l (y,v)-<S>- 1 (y,0))\ 5 ~ 1 dydv 
JGxB £ 

<kJ c k o $-l(y, 0) l^- 1 ^, v) - $- 1 (y,0)|| (d - fc)(5 ~ 1) dydv 
JGxB E 

<k x k 2 [ c k o^- 1 (y,0)\\v\\ ( - d - k) ^- 1) dydv 
J GxB e 

= k iK2 [ c k o^- 1 (y,0)dy [ lluH^K'- 1 )^, 



where Cfc(-) is as in Lemma [3l The integral over G is finite by (uniform) continuity of c k 
and 3> _1 , so we get 

m{U £ )<C [ WvW^-^-^dv, e<e Al, 
Jb £ 

for some constant C > that does not depend on e. Since B £ = eB\, the right side 
equals 

C f \\ev\\^- k ^- 1 h n2 dv = e^ d - k ^ 1 - 5 ^C [ M^W 1 ^ dv . 

JB 1 JBi 

Since the integral is over ^-dimensional space, the right side is finite provided 

(d- k){l- 5) <n 2 -2. (9) 
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But ri2 = (d — k) 2 , so Q is equivalent to [d — k)(d — k — 1 + 8) > 2, which holds for all 
k < d — 2 since <5 > 0. We conclude that there is a constant C > 0, independent of e, 
such that 

\m(U e ) < C ' £ -2-(d-k)(l-S)+n 2 _ 

e z 

Since, as we just saw, Q holds, this quantity tends to zero as e tends to zero. This finishes 
the proof of the assertion for 5 > 0. 

If 5 > 2, then Q holds also for k = d — 1, which takes care of this case as well. The only 
case that remains to prove is 5 = 2 and k = d — 1. The argument is exactly analogous to 
the one above, except that g £ is now given by 

g £ (x) = 4> £ (\7r(x)\), 

where (ft £ is the function from Lemma [H Note that = (d — k) 2 = 1, so that 7r(x) is a 
real number. The function g e is again given by (|7|). Next, instead of ([8]) we need a more 
precise estimate. Specifically, we have the inequality 

||V 5e 0r)|| <C|^(|vro$(x)|)|, 

where as before, C = sup xg [/ e ||<E>'(x)||. Hence, by the same calculations as above involving 
Lemma El we get 

/ \\Vg e (x)\\ 2 m(dx)<C 2 [ \(f>' £ (\Tro^(x)\)\ 2 \detx\dx 

= C 2 [ l^d^Dj'ldet^-^y,^!^,^®^ 

< K1 K 2 C 2 J^ca-xoQ-^y^dy \4>' £ {\v\)\ 2 \v\dv. 

By the property © of <j> e given in Lemma (U the right side tends to zero as e I 0. This 
concludes the proof. ■ 

Equipped with Theorems [9] and [TU1 we can make further progress toward proving Theo- 
rem[5J Let us start by briefly reviewing some properties of the space W 1,2 {VL, m). For any 
open set with compact closure contained in fi, we have C _1 < (detx) 5-1 < C for some 
constant C > 1 and all x € fi'. Hence 

II • \\w^{Q.',m) and II • \\w^ 2 {Q,',dx) are equivalent, (10) 

and the unweighted space W 1,2 (Q', dx) coincides with W 1,2 (Q,', m). One consequence is 
that any / £ W 1,2 (Q,m) supported on a compact set K C Q can be approximated by 
functions from C^°(Cl) via mollification: take Q' bounded with K C W C Q' C Q, so that 
/ € W 1,2 (Q', dx), and use standard mollification results for this space. 

Another useful consequence of (]10p is that W 1,2 ($l,m) is stable under truncation: If 
/ G VF 1,2 (0, m) then g = f Al has a weak derivative with = Z)/l|j <1 } (in which case, 
of course, ||5 f ||iy 1 . 2 (n,m) ^ Il/llvi /1 ' 2 (r2,m))- To see this, pick any test function ip G C£°(£l) 
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with support K C f2. Now choose Q,' as above, and apply known truncation results for 
W^ 2 {p!,dx) to get 

(f(x)M)ip(x)dx = - / (Df(x)l {f{x)<1} )»Vip(x)dx. 

Since (p vanishes outside f2', the above relation holds with replacing £1' . This proves the 
claim. It is then easy to show that 

the bounded elements of W 1,2 (£l,m) are dense. (11) 

After this digression we continue with two technical lemmas. The first shows that any 
element of W 1,2 {£l,m) is close to some element whose support is bounded away from 
Ufc<d-2-^fc (and from dQ = Uk<d-iMk in case 5 > 2.) 

Lemma 5 Let 5 > and pick f E W ' (Q,m). For any e > there is g E W 1,2 (Cl,m) 
with bounded support such that \\f — g\\w 1 ' 2 {Q,m) < £ an d g = on U HQ for some open 
set U C M d containing If 5 >2, U can be chosen to contain dQ,. 

Proof. It is enough to consider positive / with bounded support, and by (lllh we may 
assume / is bounded. By Theorem [TO] there are open neighborhoods U n of ^k<d-2^k 
(respectively of <9il if S > 2) and elements h n E W 1,2 (M. d , m) such that h n = 1 on U n 
and h n — > in M^ 1,2 (r2, m). Replacing /i n by (h n A 1) V if necessary, we may assume 
< hn < 1. Now define 

/« = (1 - K)f. 

Then 

ll/n " /llwi.»(n,m) = J [\K{x)f{x)\ 2 + \\h n {x)Df{x) + f{x)Dh n {x)\\ 2 ) m(dx) 
<2||/||Lll^|| 2 wl , 2(c , m) + 2^|/ ln (x)| 2 ||^/(x)||V(^)- 

Fix any e > 0. Let A K = {x £ E : \\Df(x)\\ < k}, and pick k large enough that 
f Ac \\Df{x)\\ 2 m(dx) < e/2. With this k, 

h n (x)\ 2 \\Df(x)\\ 2 m(dx)< [ \h n (x)\ 2 \\Df(x)\\ 2 m(dx)+ [ \\Df(x)\\ 2 m(dx) 
n Ja k JAi 



< k I \h n (x)\ z m(dx) + e/2. 
Jn 

It follows that for n sufficiently large, 

ll/n _ /llw 1 .2(n,m) < £ - 

We may now take g = f n - ■ 

The proof of the next lemma uses a standard partition of unity argument (mimicking the 
proof of [H Theorem 3.22]) together with a powerful extension theorem due to Chua [5] 
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for weighted Sobolev spaces with Muckenhoupt weights. In particular, therefore, we will 
rely on the Muckenhoupt property of Idetxj 5 " 1 for < 5 < 2. One delicate point is to 
choose the covering sets so that Chua's theorem becomes applicable. Specifically, their 
intersections with fi must be so-called (e, 5) -domains. Here the tube segments discussed 
previously will again be useful. 

Lemma 6 Let < 5 < 2 and define fi = fi \ {x G M d : ranks; < d - 2}. Then 

C c °°(fi) is dense in PV" 1>2 (fi,m). 

Proof. By Lemma [5] it suffices to approximate elements / G W 1,2 (£l, m) with bounded 
support and with / = on U Pi fi for an open subset U C M d containing Ufc<d_2-/Vffc. 
Define 

K = closure in M d of {x G fi : / 0}, 

and note that X is compact. Also, by shrinking U if necessary we may assume that 
[/ fl A - = t). Now, for every x G M^-i, let C/ x be an open tube segment around x, thin 
enough so that rank y > d — 1 for all y £ U x . Note that the set 

F = K\U xeMd _ 1 U x 

is compact in M" 1 and is contained in fi. Hence there is an open set Uq such that F C 
Uq C 0. Next, since K is compact, we may choose finitely many of the U x , say U±, . . . , U n , 
such that K C U" =0 ?7j. There are also open sets Vb, . . . , V n with 7, C (7j for each i and 
such that K C U™ =0 Vi. Let {ipi : i = 0, . . . , n} be a smooth partition of unity subordinate 
to Vb, . . . , V n , and define fi = ipif. Fix any e > and suppose we can find fa G C£°(fi) 
such that || /j — 4>i\w^ 2 (Ci,m) < + n )- Then, letting = Yli&i we nave S C£°(fi) 
with 

n 

as desired. We thus focus on finding suitable (pi. For i = one can simply take a 
mollification of /o, thanks to (I10p . For i G {l,...,n}, we define fij = n f/j. Then 
/i G W 1,2 (fij,m), and it is readily verified that fij is a Lipschitz domain. Indeed, by 
construction it is the image of a cylinder — the product of a Euclidean ball and a (short) 
open line segment — under a bi-Lipschitz map. Since | det is a Muckenhoupt A2 weight 
by Theorem [9] (and since every Lipschitz domain is an (e, <5)-domain, see [15} p. 73]), the 
extension theorem of Chua [3 Theorem 1.1] yields an element gi G W 1,2 (WL d ,m) with 
9% = fi on Since fi = outside V$, we may multiply gi by a smooth cutoff function 
that is one on Vi and zero on Uf, and thus assume that gi = fi on all of fi. The desired 
4>i G C^°(fi) is now obtained by mollification of gi. Indeed, the support of (pi can be made 
bounded away from L)k<d-2Mk, an d we have 

\\<j>i - fi\\w 1 ' 2 (E,m) = \\4>i ~ 9i\\w 1 ' 2 (E,m.) < W^i ~ 9i II W 1 ' 2 (M d ,m) > 

where the right side can be made arbitrarily small. ■ 
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We are now finally ready to prove our main result leading to Markov uniqueness of the 
BESM semigroup. 

Proof of Theorem [8]. In view of ([5]), it is clear that in order to prove (i) (re- 
spectively (ii)) it suffices to prove W l ^(Q,,m) C H2\n,m) (respectively W^fi.m) C 
H^ 2 (n,m).) 

We start with part (ii), and pick / G W 1,2 (Q,,m). By Lemma / can be approximated 
by some g G W 1,2 (Q,m) whose support is bounded and bounded away from dQ. By 
mollification (relying on (110p ). g can in turn be approximated by h G C£°(f2). Since 
HQ ,2 (Q,m) is the completion of the set of such functions, we obtain / G HQ ,2 (Q,m), as 
desired. 

Now consider part (i). By Lemma [6] the result follows if any given / G C^°(£l) can be 
approximated in W 1,2 (Q, m) by g G C 2 j_. For this we can apply the same modified version 
of [HJ Theorem 5.1] as in the proof of Proposition [2] above. However, since / is zero in a 
neighborhood of L)k<d-2Mk, there is a simple direct argument that we now give in order 
to make the presentation self-contained. 

For each e > 0, let 4> £ G C 2 (IR-|_) satisfy the following properties: 
(i) £ (i) = for t > e, 
(ii) iz^W and M) 

are bounded in t (where the bound may depend on e), 
(iii) \<p e (t)\ < e and \4>' e it)\ < 2 for a11 * > °- 

This is clearly possible when e = 1, and hence also in general by setting 4> e (t) = £<j)\(t/e). 
Let K denote the support of /. With g(x) = ^vd'T^p^ we ^ ave 

Vf(x) = g(x)Vdet(x) + G(x), (12) 

where G(x) • Vdet(x) = for all x. Note that Vdet(x) = adj x T ^ for all x G K, 
so that (7 is well-defined if we set it to zero outside K. In fact, g G C£°(E). For e > 0, 
consider the function 

h £ (x) = f(x) - (p £ (detx)g(x). 

Straightforward calculations using the chain and product rules together with (fT2j) yield 
the following two expressions for V/i £ : 

Vh £ (x) = Vf(x) - (f>' e (detx)g(x)Vdet(x) - <p £ (detx)Vg(x) 

= (l-^(detx)) Vf{x) + (/)' £ (detx)G(x) - <j> e (debx)Vg(x). 

We first show that h £ approximates / in VF 1,2 (i7,m) if e is chosen sufficiently small. But 
this is an immediate consequence of the inequalities 

\f{x)-h £ {x)\<e\g(x)\ 

and 

||V/(s)-VMx)|| < (2|5(x)|Vdet(x)|| + ||V 5 (x)||)l [0ie] (detx), 
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together with the fact that m{{x G K : < det < e}) tends to zero as e 4- 0. It thus 
suffices to check that h £ G C%±- Clearly h £ G C^(Q). Moreover, the second of the above 
expressions for V/i £ together with the orthogonality of Vdet(x) and G(x), as well as the 
fact that x~ T = Vdet(x)/ det (a;), yield 

x~ T • Vh £ {x) = 1 ~^} de \ x) Vdet(x) . V/(x) - ^ d6tx) V det(x) • Vg(x). 
(p £ (detx) detx 

The properties of <p £ imply that the right side is bounded, as required. Part (ii) of the 
theorem is thus proved. 

It remains to prove Markov uniqueness. But this follows directly from the basic criterion 
for Markov uniqueness given in |9l Corollary 3.2], since for any 5 > we have W 1,2 (Q., m) = 
i^' 2 (Sl,m). ■ 

Remark 11 Using results in [8], the space C^ ± is dense in C^{E) with respect to the 
norm || • ||w /1 ' 2 (r2,m)- A n alternative approach to proving Theorem[8] (i) would therefore be 
to show directly that C^(E) is dense in W 1,2 {£1, m), for example by showing that is an 
(e, 5)-domain and then apply Chua's extension theorem. Proving the (e, 5) property does 
not appear to be straightforward — one obstruction is that £1 does not lie on one side of its 
boundary, as discussed in the Introduction. 



5 Proof of the integration by parts formula 

In this section we give a proof of the integration by parts formula, Theorem [2j which we 
now restate for the reader's convenience: 

Assume that 5 > 0, and let f £ Cl(E) and G £ C X {E; M d ). If 5 < 1, assume that G(x) is 
tangent to dE at x for all x G dE. If 5 < 1, assume in addition that G(x) • x~ T is locally 
bounded. Then 

(Vf,G) = (f,-V*G). (13) 



Throughout the proof, let K be the (compact) support of /. For e > 0, define 

Vdet(x) 

U £ = {x G M d : det x > e} and v(x) = -—— — x G U - 

For e > 0, U £ has smooth boundary with outward unit normal v{x) at x G dll £ . For 
any smooth function h : Uq — > R such that the integrals are well-defined, the standard 
integration by parts formula yields, for each e > 0, 



V/(x) • G(x)h(x)dx = [ f{x)h{x)G(x)»v{x)da £ {x) 

- [ f(x)V.(Gh)(x)dx, (14) 
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where a e denotes the surface area measure on dU e . 

Case 1: 6 > 1. Take h(x) = det(x) 5 - 1 . As e | 0, the left side of flH) tends to J E Vf(x) • 
G(x)m(dx) by dominated convergence. Let C > be such that \ f(x)\ \\G(x)\\ < C for all 
x £ K. The absolute value of the boundary term is then dominated by 

Ce 5 - 1 f da £ {x), 
JdU £ nK 

using also that h(x) = e* 5-1 for x G dU £ . It is easy to see that cr e (K) remains bounded 
as e 4 0, so we conclude that the boundary term vanishes in the limit. Consider now the 
second term on the right side of (|14p . The product rule yields 

V»(G7t) = hV»G + G» V/i. 

By dominated convergence, fjj f(x)VG(x)h(x)dx — > f E f(x)VG(x)m(dx). Moreover, 
G(x)*Vh(x) = (6— l)G(x)» V det(x) det(x) 5-2 , and since det(x) s ~ 2 dx is a Radon measure 
due to Theorem Q] and the fact that 5 > 1 , we may again use dominated convergence to 
get 

/ f(x)G(x) • Vh(x)dx -> (S - 1) / f(x)G(x)»x~ T m(dx). 

Ju e JE 

(Here we used the equality V det(x) det(x) 5 ~ 2 dx = x~ T m(dx).) Assembling the pieces 
gives the desired formula (|13p . 

Case 2: 6 = 1. We again take h(x) = det(x) s ~ 1 = 1. Except for the boundary term, 
everything works as in the case 6 > 1, if we just note that V/i = 0. Letting C be a bound 
on |/(x)| over K, the boundary term is dominated by 

Ca £ {K) sup G{x)*v{x). 

xeu e nK 

Using that G(x) is tangent to dE at every x G it is not hard to show that the supremum 
tends to zero. Hence (fT3|) is established. 

Case 3: 6 < 1. Things are now a bit more complicated due to the fact that det(x) 5 " 1 
blows up at dE. To get around this, for each n let r n be a smooth, nondecreasing function 
satisfying the following properties: 

T n (t) < t A n, r n (t) = i for t < n — 1, r n {t) = n for t > n + 1 

< < 1, r n (i) t * as n -)• oo. 

In (JTIJ) we now take /i = h n , where h n = r n ow and w(x) = det(x) s ~ 1 . We first hold n fixed 
and let e ! 0. The left side of (fTl|) converges to f E S/f(x) • G(x)h n (x)dx by dominated 
convergence. The boundary term on the right side will vanish by the same argument as 
in the case 6 = 1. The integrand in the second term on the right side is in fact bounded, 
since by the properties of r n , 

X7h n (x) = (6- l)«ou,(x))Vdet(x)det(x) 5 - 2 < (6 - 1)V det(x)(n + 1)^. 
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Dominated convergence gives the limit J E f(x)V • (Gh n )(x)dx. Combining these results 
and applying the product rule gives the formula 

Vf(x)»G(x)h n (x)dx = - j f{x)(v •G{x)h n {x) + G{x)»Vh n {x)\dx. 

J E 

The final step is to send n to infinity. The left side converges to f E V/(x) • G(x)m(dx) 
by dominated convergence, since h n = r n o h f h. Similarly, J E f(x)V • G{x)h n {x){dx) — > 
L /(x)V • G{x)m{dx). Finally, 



A' 



= (5-1) 






'if 




<(8-l) 


L 


G(x) • x 








' • x~ T is 


bounded on 



G(x) 



Vdet(x) 



det(x) 



m(dx) 



m(dx) 



convergence yields J E f(x)G(x) • Vh n (x)(dx) — > f E f(x)G(x) • Vh(x)(dx), and hence the 
result. 



6 Proof of the Muckenhoupt A\ property 

This section is devoted to the proof of Part (i) of Theorem [9j That is, we fix a G (—1,0] 
and define 

w(x) = | det x\ a , 

where we set w(x) = +oo if detx = 0. We will prove the existence of a constant C > 0, 
depending only on d and a, such that the inequality 

— — / w(x)dx < C inf wix) (15) 

\B\ J B x€B 

holds for every ball B C M d . Define the set D^. of diagonal matrices with nonnegative 
and ordered diagonal elements, 

= {Diag(o-) : a G R rf , a x > a 2 > ■ ■ ■ > a d > 0}. 

The open ball centered at x £ M d with radius r > is denoted by B(x, r). Its intersection 
with the nonsingular matrices is denoted by B*(x,r). That is, 

B(x, r) = {y £ M d : \\x - y\\ < r}, B*(x, r) = {y G B(x, r) : det y / 0}. 

The proof of the Muckenhoupt property is quite long and involved (but nonetheless mostly 
elementary), due to the relatively complicated geometrical structure of the set dE, which 
is where the weight function becomes singular. The main idea is to change variables and 
integrate over the product space 0(d) x T(d) instead of M d . Unfortunately, balls in Nl d 
do not always map to balls (or comparable shapes) in 0(d) x T(d), and this is where the 
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main complications arise. The resolution to this issue resides in Lemma [8] below, which 
relies on a detailed analysis of the mapping taking x to its Q-R-decomposition. 



We start with a lemma that establishes an inequality similar to (115f) , where the balls B are 
replaced by sets of the form U -K, with U C 0(d) measurable and K C T(d) a cube. 

Lemma 7 Let — 1 < a < 0. Then there is a constant C\ > 0, depending only on d and 
a, such that the inequality 

/ w(x)dx < C\\U ■ K\ inf w(x) 
Ju-K 'xeU-K 

holds for any measurable subset U C 0(d) and any cube K C T(d). 

Proof. Pick a cube K = {R € T(d) : Rij G Iij, i < j}, where the Iij are intervals, and 
let U C 0(d) be measurable. By Lemma [U we have 

f w(x)dx = n(U) f J] Ri~ i+a dR = fi(U) J] Ihj] H f t d ~ l+a dt, 

JU-K JK i=1 i<: - i=1 Jlu 

and similarly \U ■ K\ = n(U) \hj\\[* =1 J T .. t^dt. We also have 

d d 
inf w(x) = inf TlRfi = TT inf t a . 

X&U-K RG/C 11 1L t£lii 

Hence the result follows from the following Claim: 

Let a G (—1,0] and f3 > 0. Then there is a constant C a ^ such that for every interval 
I C (0, oo), we have 

J^t a+ Pdt<C a ^ Jt p dt inft a . 
To prove the Claim it suffices to consider / = (a, b) with b > a > 0. We obtain: 

Ji a + /3 + lV Kb) ) 

< (bP +1 -aP +1 ) b a 

~ a + /3 + 1 V J 

= ftPM inft", 

a + + 1 7/ te/ 

as required. ■ 

Consider now balls B(S, r), where the diagonal elements of S G are either "large" 
(comparable to the radius r) or zero. The following result reduces the proof that (|15p holds 
for balls of this form to an application of Lemma [7J In the statement of condition (|16j) 
below, we use the convention that uq = oo and that i runs over {0, . . . , d}. 
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Lemma 8 Suppose E G and r > satisfy the following property, where a G R rf is £/ie 
vector of diagonal elements of £ : 

There is an index n G {0, 1, . . . , d} such that . . 

o~i > (5 + 18d)r /or a// i < n, and o~i = /or a// i > n. 

Then there is a measurable subset U C 0(d) and a cube K C T(d) such that 

B^,r) CU ■ K C B^,C 2 r), 

where C 2 is a positive constant that only depends on d. 

Proof. The problem of finding the advertised constant C 2 can be reduced to proving the 
following Claim, where e\, . . . , denote the canonical unit vectors in ~R d : 

There is a constant O3, depending only on d, such that the following holds: For any 
x G £?*(E,r), let x = QR be its QR- decomposition, and let qi,...,qd be the columns of Q. 
Then the inequalities \\R — S|| < C%r and \qi — | < ra^ x Cs hold for all % G {1, . . . ,n}, 
where n is the index from condition (j 16f) . 



Let us show how the statement of the lemma follows from this claim. Define K to be the 
cube in T(d) centered at £ with side 2C^r, i.e. 

K = {Re T(d) : max |^ - Sy| < 3 r}, 

and let U C 0(d) be given by 

f r 

*7 = < Q = (qi, ■ ■ ■ ,qd) € O(d) : |<7i - e» < — 3 ,i = 1, . . . ,n 

The Claim then directly implies .B*(£,r) C U ■ K. Next, for any x = QR £ U ■ K the 
triangle inequality and the definition of K yield 



\\x - E|| < - E|| + 11(0 - I)E|| < y/d(d + l)/20 2 r + ||(Q - I)E||. 

Furthermore, since o"j = for i > n, we have ||(Q — I)E|| 2 = a\\q\ — e\\ 2 + . . . <r^\q n — e n \ 2 . 
We then deduce from the Claim that \\(Q — I)E|| < ^/nC-jr, and consequently 

Hz - E|| < C 2 r, where 2 = (^d(d + l)/2 + O3. 



We are thus left with proving the Claim. Since it is vacuously true for n = 0, we can assume 
n > 1. The proof relies on a rather careful analysis of the Gram-Schmidt orthogonalization 
procedure for obtaining the 0-^-decomposition °f a generic matrix x G _B*(E,r), so we 
briefly recall this procedure. To improve readability, we temporarily (for this proof only) 
adopt the notation (y, z) = y T z for y, z G H d . Fix x G B*(£, r) and let x\, . . . ,Xd be the 
columns of x. To obtain the Q/2-decomposition of x, one defines 

Ul 

ui = xi, qi = ] — r, 
Pi 
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and, if gi, . . . , <fy-_i have been defined, 

3-1 



Uj — Xj 



^2(qi,Xj)qi, qj = ^-. (17) 



i=l 



The vectors qi,...,qd obtained in this way are the columns of Q, and R is given by 
Rij = (qi,Xj), i < j. 

We now proceed with the proof of the Claim. Recall that e±, . . . , e<i are the canonical unit 
vectors in R rf . Since x £ B(E, r), we have Xj = cr^ei + hi, where hi is a vector in R rf with 
\hi\ < r. Also let a = 5 + 18cZ denote the constant appearing in condition (I16|) . 



Fix j £ {1, . . . , n}, and suppose we have proved the following: 

For all i < J — 1 and all A; > i, \R ik \ < 3r. (18) 
Then (|17p and the inequalities \xj\ > cr,- — |/t,-| > a a — r imply 

3-1 

\uj\ ><Tj-r-Y, \Rij\ > 0j ~ r(l + 3(j - 1)). (19) 

i=l 

Moreover, for k > j we use (fTKj) and the fact that \(xj, Xk)\ < o^r + cr^r + r 2 to get 

i-i 

|(«i,a:fc>| < K^-^fe)! +^1^-1 \Ri k \ < a jr + a k r + r 2 {l + 9(j - 1)). 

i=i 



It follows that 



_ |(Mj,x fc )| gj + a k + r(l + 9(j - 1)) 
' K-l " ^-r(l + 3(i-l)) 

l + fffc/ffj + Cr/ffj-Kl + SO'-l)) 



< r 



l-(r/^)(l + 3(j-l)) 

2 + q- 1 (l + 9(i-l)) 
l-a-i(l + 3(i-l))' 



where in the last step we used that a k < 0j (since fc > j) and > ar (since j < n). 
Since a = 5 + 18d, the right side is smaller than 3r, as one readily verifies. We deduce 
that (USD holds with j replaced by j + 1, and since it is vacuously true for j = 1 it follows 
by induction that it holds for all j £ {1, . . . , n + 1}. 

We now use this result to bound \Rjj — Oj\ for j £ {1, . . . , n}. To this end, write 

\ x j\ 2 ~ a j\ u j\ \ + %3 ~ l) 7 " 2 



\ R 33 ~ a j\ ~ J— 



3-1 

|2 



i=i 

using that (qi,Xj) 2 = \Rij\ 2 < 9r 2 due to (fTB~l) . Moreover, we have 

= \a 2 — <Tj\itj\ + \hj\ 2 + 2aj(ej, hj)\ < Oj \oj — \uj\\ + r 2 + 2ajr, 



< 
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and by the reverse triangle inequality, 

\&j — \ < \&jej — Uj\ < r + 3(j — l)r. 
Assembling the pieces and using the bound (jTUJ) gives 



(20) 



\Rjj -<tj\ < 



a jr {l + 3{j - 1)) + r 2 + 2a jr + 9(j - 1> 
aj - r(l + 3(j - 1)) 



Dividing the numerator and denominator by <x,- and using that Oj > ar, we finally arrive 
at 



\Rjj-°j\ < r 



< r- 



1 + 3(j - 1) + a" 1 + 2 + q-^C? - 1) 

l-a-^l + SCi-l)) 
S + SO'-^ + q-^l + oq-l)) 

l-a-i(l + 3(i-l)) 
14 + 12<i 



The only elements of i? that remain to analyze are Rij for i < j and j > n. But = /ij 
for these j, so = |(oj, Xj)| < \hj\ < r. We are now able to estimate \\R — S|| as 

follows: 



\R- Ell 



l^'l 



< r 



14+12d , ra(n-l) 

x d + 3 x ^ + d(d - n) 



3 2 

A bound solely in terms of d is then easily obtained. For instance, we may take 

,"14 + 12d (d-1) 
|i? - S < rd ( + 3 x v n ' + d 



(21) 



Let us now focus on bounding \qj — eA, j £ {1, . . . , n}. The calculations are similar to the 
ones used to bound \Rjj — Uj\ above, but slightly simpler. We have 



m ~ e j\ 



j'-i 



t=l 



< 



(T 



>• - Kl| + r + 3(i - 1)) 



< (r + 3(j-l)r , 
l«il v ' 

using (j20|) in the last step. Using again (|19p together with <jj > ar, 

r 2 + 6(7-1) r 5 + 15d 

o,- — e,- < — x — ; — < — x . 

m Jl ~ Uj 1 - q-!(l + 3(j - 1)) Oj 2 

The Claim, and hence the lemma, is now proved, if for C3 we take the maximum of 
2(5 + 15d) and the constant in (]2ip . ■ 



The next result implies that the proof of (|15j) for any ball with diagonal center reduces to 
an application of Lemma [H 
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Lemma 9 Pick any E = Diag(<r) G and r > 0. There is a matrix E' = Diag(o"') G 
and a real number r' > t/iai satisfy the condition (I16p (with a replaced by a' , and r by 
r' ), such that 

5(E, r) C £(E', r ) and r' < (6 + 18d) d r. 

Proof. Let a > be a constant to be determined later. Suppose for some index i G 
{0, . . . , d— 1}, we have Od~i > a(l+a)V. Let k be the smallest such index, and define 

a' = (<7i,...,<j d _ fc ,0,...,0), r' = r(l + a) fc . 

Then, since cr^-i < a(l + a)V for all i < k, 

\a - a'\ < o-d-k+i H K cr rf 

< r fa(l + a) k ~ l H h a(l + a) + a 

= r(l + a) fc - r 

By the triangle inequality, B(E, r) C B(T,', r'), where E' = Diag(cr'). Setting a = 5 + 18d, 
we see that E', r' satisfy condition (|16p with n = d — k. 



It remains to consider the case where < a(l + a)V for all i G {0, . . . , d — 1}. Then 
any x G 2?(E,r) satisfies 

IHI <r + r (a(l + a)^" 1 H h a(l + a) + a) = r(l + a) d , 



so B(E,r) C 5(0, r(l + a) rf )- With r' = r(l + a) rf and n = 0, condition (fT6j) is again 
satisfied for a = 5 + 18d. This finishes the proof. ■ 

The proof of (|15p is now straightforward. Indeed, pick any ball B = B(x,r), and let 
x = UZV T be a singular value decomposition of x. Then B(x,r) = U ■ £?(E, r) • y T , and 
together with the invariance of Lebesgue measure under orthogonal transformations and 
the fact that dety = det(U T yV) for any y G M d , this leads to the equalities 



w(x)dx = / w(x)dx, 

B(x,r) JB(E,r) 

\B(x,r)\ = |£(E,r)|, 



inf w(x) = inf w(x). 
x£B(x,r) xeB(E,r) 



Consequently (and using that dE is a nullset), it suffices to prove f)15|) for i? replaced by 
-B* = -B*(E,r) with E G Di. We then have the following chain of inequalities, where we 
set Bl = 2?*(E',r') with E' and r' from Lemma El and where U, K, and C2 are obtained 
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(B.cBicU- K) 
(Lemma [7]) 
(B*cU- K) 
(U ■ K c £*(£' \C 2 r')) 
(r' < (6 + \M) d r) 
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